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■ Recent striking lattice results on strong interaction and bound 

, states above T c can be explained by the nonperturbative QQ poten- 

| tial, predicted more than a decade ago in the framework of the field 

correlator method. Explicit expressions and quantitative estimates are 
given for Polyakov loop correlators in comparison with lattice data. 
New theoretical predictions for glueballs and baryons above T c are 



qj\ given. 



1 Introduction 

X 

There is a growing understanding nowadays that nonperturbative dynamics 
plays important role in the deconfinement phase, for reviews and references 
sec [1. 

An additional part of this understanding, not contained in pQ, is the re- 
alization of the fact, that at T c < T < 2T C , the colormagnetic fields are 
as strong as they are in the confinement phase, (where colormagnetic and 
colorelectric fields are of the same order) and become even stronger above 
2T C . More than a decade ago the author has argued [21 E] that the deconfine- 
ment phase transition is the transition from the color-electric confinement 
phase to the colormagnetic phase, confining in 3d. This observation was sup- 
ported theoretically by the calculation of T c Ej and on the lattice by the 
calculation of the spacial string tension at T > T c 
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In 1991 the author has found [5] that colorelectric fields also survive the 
deconfinement transition in the form of potential V\(r), which can support 
QQ bound states in some temperature interval T c < T < To, while quarks 
acquire self-energy parts equal to |Vi(oo). As will be shown below this pre- 
dicted picture is fully supported by recent numerous lattice experiments jH]- 
[T3| (for a review see (101), where QQ bound states have been discovered. At 
the same time the evidence for Vi(r) and selfenergies ^Vi(oo) has been also 
obtained on the lattice in the form of Polyakov loop averages and of free and 
internal energies above T c [TUl HTJ U2] • The light quark (m q ~ m s )qq bound 
states have also been observed in [T3j . 

The theory used in jH] and below is based on the powerful Method of 
Field Correlators MFC [T3|, (for a review see where the basic dynamic 
ingredients are the field correlators (trF A11 j, 1 (xi)...F AtnI/n (a; n )) 1 . It was shown 
later [TH] that the lowest quadratic (so-called Gaussian) correlator explains 
more than 90% of all dynamics and it will be considered in what follows. 
The quadratric correlator consists of two scalar form-factors, D and D\\ 

D^x, 0) = ^tr(F flv (x)^(x,0)F Xa (0)) = 
= {SpxHvo - 5 lxa 5 vX )D{x) + \ -^{xxSua - x a 5 uX ) + (//A ucr) D x {x) 

(1) 

which produce the following static potential [TZJ between heavy quarks at 
zero temperature (obtained from the Wilson loop r x t with t — > oo) 



V(r) = 2r J r dX J™ drD{ VW+r ^) + J r XdX / °° rfr[-2D(VF+^)+ 

+D 1 (VX T T^)] = V D (r) + V 1 (r). { ] 

One can notice that linear confinement part of potential, Vd = crR, is 
due to correlator D(x), a = | J D(x)d 2 x. 

At T > one should distinguish between electric and magnetic corre- 
lators, D E (x), Df (x), and D H (x), D^(x) and correspondingly between a^ E > 
and a^ H \ It was argued in [21 E] that the principle of minimality of free 
energy requires D E and electric confinement , a^ E \ to vanish, while color- 
magnetic correlators, D H (x),D 1 H (x) should stay roughly unchanged at least 
up to 2T C . 

Several years later in detailed studies on the lattice in [TH] these state- 
ments have been confirmed, and indeed magnetic correlators do not change 
at 1.5T C > T > T c while D E (x) vanishes in vicinity of T c . 

1- we omit for simplicity parallel transporters necessary to maintain gauge 

invariance, see ^H] for details. 
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Not much was said about the second electric correlator Df(x), since in 
the parametrization of 18 it was found to be smaller than D E (x) and hence 
not so important at T < T c . 

Meanwhile a lot of information was being accumulated on the lattice. 
First of all, the Polyakov loop averages already imply the presence of strong 
electric fields above T c , and the main point is that those can not be reduced 
to the perturbative electric and magnetic, (see the analysis in [TT| IT2]). 

Recently a detailed analysis of Polyakov loop correlators was done by 
the Bielefeld group an d the singlet free energy Fi(r,T) and internal 

energy Ui(r,T) were calculated at T < T c and at T > T c . In the latter case 
.Fi(r, T) was found to saturate at large r at the values of the order of several 
hundred MeV (e.g. for T = 1.2T C the value of F 1 (oo,T) found in |T2] is 
around 0.7^/o : , while the internal energy is around 3T C ) and this fact cannot 
be explained by perturbative contributions alone - we consider it as the most 
striking revelation of nonperturbative electric fields above T c . 

At the same time several groups have calculated the so-called spectral 
function of heavy [^-[10] and light quarkonia [T3j at T > T c . In all cases 
sharp peaks have been observed, corresponding to the ground state levels of 
cc at L = 0, 1 and of light (m q « m s ) quarkonia in V, A, S, PS channels. In 
both heavy and light cases the peaks are possibly displaced as compared to 
T = positions and apparently almost degenerate in different nn channels. 

All these facts cannot be explained in the framework of the commonly 
accepted perturbative quark-gluon plasma and call for a new understanding 
of the nonperturbative physics at T > T c . In what follows we shall argue 
following |5] that at T > T c not only nonperturbative magnetic fields, but also 
strong nonperturbative electric fields are present, which can be calculated in 
MFC and explain the observed data. 

2 Dynamics of Polyakov loops and the corre- 
lator D\ 

In this section we consider the Polyakov loop and apply to it the nonabelian 
Stokes theorem and Gaussian approximation, taking first the loop as a circle 
on the plane and making limiting process with the cone surface S inside loop 
and finally transforming cone into the cylinder by tending the vertex of the 
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cone to infinity 2 . In doing so we are writing the nonabelian Stokes theorem 
and cluster expansion for the surface S which is transformed from the cone to 
the (half) cylinder surface. As a result one has for the Polyakov loop average 



(3) 



L = j^trP exp (-| Jgda^u) j s d<T p \(v)D mp x{u,v) 
= exp \-\ J 1/T drti /T dr> / °° ^D^ 2 + (r - r') 2 )} . 



In obtaining we have omitted the contribution of D(x) in D^ vp x, since 
this would cause vanishing of L in the limiting process described above due 
to the infinite cone surface S. This exactly corresponds to vanishing of L 
in the confinement region, observed on the lattice. Therefore the result (jSJ) 
refers to the deconfinement phase, T > T c . 

As it is known from lattice ^H] and analytic calculations, -Di(x) [T§] 
exponentially falls off at large x as exp(— Mix), with Mi > 1 GeV and for 
T <C Mi one can approximate © as follows 3 

L = exp (-| J 1/T dv (1 - v) J °° (P^VFT^)) « 

w exp (-5^(00)), \/ 1 (r,T) = / 1/T dz/(l-^)/o r ^i(v / e+^) ' 

(4) 

We turn now to the correlator of Polyakov loops following notations from 
jllj . Using the same limiting procedure as for one Polyakov loop, one can 
apply it to the correlator (trL^trLy) = P(x — y), tr = j^tr, representing 
the loops Lx and Ly as two concentric loops on the cylinder separated by 
the distance |x — y| along its axis, the cylinder obtained in the limiting 
procedure from the cone with the vertex tending to infinity. One can apply 
in this situation the same formalism as was used in [201 for the case of the 
vacuum average of two Wilson loops. For opposite orientation of loops using 
Eqs. (21-28) from [20] one arrives at the familiar form found in 



P ( \ 1 ( P ^ T ) \ ■ N c~ l ( F M \ _, i 



c \ / c 



(5) 



2 In doing so one is changing topology of the surface and as a result loses the Z(N) 
subgroup of SU(N). This however does not influence our results as long as one is remaining 
in the j — sector of Z(N) broken vacua (see last ref. of for more discussion of Z(N)) 

3 The correlators D, D\ in J2J in principle should be taken in the periodic form, as was 
suggested in . However for T < 2T C this modification brings additional terms of the 
order of exp(— Mi/T) which are neglected below. 
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In the Appendix two different ways of derivation of Eq.(|5|) are given, with 
the result 

F 1 (r,T) = V 1 (r,T) + V D (r,T), (6) 
exp(-F 8 (r,T)/T) = L adj (T) ex V {-(V D (r,T) - ^(r, T))/T}. (7) 
Here 

L ad] = exp(-(^ D (r*) + ^(oo, T))/T) (8) 

is the vacuum average of the adjoint Polyakov loop, which vanishes in the 
leading approximation in the confinement phase, as it is explained in the 
Appendix, and nonzero when gluon loops are taken into account, in which 
case \V D (r\T) ~ Mi and L adj ~ exp(-M x /T) < 1. 

The suppression of exp(— F 8 /T) in our approach in the confinement phase 
has thus the same origin as the strong damping of the adjoint Polyakov loop 
in that phase [23] and the persistence of the Casimir scaling for adjoint 
static potential in the interval < r < 1.2 fm (see |21j for discussion and 
references) . 

It is clear that in the deconfinement phase with D = 0, Vp = one 
has only Vi(r) in both F\ and F 8 , and all these quantities are finite (after 
the renormalization of the perturbative divergencies specific for the fixed 
contours, which are discussed in section 3). Thus in the deconfined phase 
one can write 

p( x _ y ) = e -*«/T = I e -V l( r)/T + 8 e -(|v l( oo)-|y l( r))/T (g) 

9 9 

where Vi(r) and Vi(oo) are renormalized. It is clear from @ that at small 
r one has lim r ^ (Fg ? (r) — Vi(r)) — > TlnN^ as was noticed and measured in 

nn. 

At this point one should stress the difference between the genuine free 
energy Fi(r, T),i — 1.8, which is measured with some accuracy on the lattice, 
and the calculated above Fi(r,T). It is clear that Fi do not contain the 
contribution due to excitation of QQ and gluon degrees of freedom existing 
at finite T. The latter is contained in the free energy i<\(r, T) and in the 
internal energy, which we denote £/»(r, T) = Fi + ST to distinguish from our 
Vi(r,T), since they are not equal. 

In general for nonzero temperature and comparing to the lattice data 
on heavy-quark potential one should have in mind, that temperature effects 
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might be of two kinds. First, the intrinsic temperature dependence due to 
changing of the vacuum structure and the vacuum correlators and hence of 
our potentials F(r,T). Second, the physical quantities like Fi(r,T),Ui(r,T) 
are thermal averages over all excited states, e.g. 



e 



-fi(r,T)/T _ V- p -E n {r,T)/T 



£e"W (10) 



U 1 {r,T)=Y J E n e~ E ^T)IT (11) 

n 

One can associate F\(r,T) = E (r,T), while the structure of excited spec- 
trum can be traced in the temperature dependence of F\ and U\. E.g. as- 
suming in the confinement phase the string-like spectrum and multiplicity for 
the multihybrid spectrum with two static quarks, E n = or + im/r,n = 1, 2, .. 
and multiplicity p(m) = exp(m/m )6(m — mi),m = im/r, one arrives at 

F 1 (r,T) = ar + m 1 (l-T/m )-Tln(--—) (12) 



T m 

U 1 (r,T) = ar + m 1 + T/(l-T/m ) (13) 

The increase of Ui(r,T) below T c in the quenched case was indeed observed 
in lattice calculations (see Fig. 3 of [T^]). 

Above T c one can see in lattice data ^0] the striking drop of entropy 
5*1(00, T) and Ui(oo,T) in the region T c < T < 1.2T C which can be possibly 
explained again by the multihybrid states occurring due to the potential 
Vi(r, T) = Fi(r,T) connecting quarks and gluons, and assuming that the 
magnitude of Vx(r,T) decreases with temperature passing at T ~ 1.05T C the 
critical value enabling to bind those states of high multiplicity. In this way 
one assumes that both below and above T c in quenched and unquenched cases 
the dominant (in entropy) configuration is the gluon chain connecting Q and 
Q with gluons bound together by confining string (below T c ) and potential 
V\ (above T c ). Thus the comparison to the lattice data on Fi,U% needs the 
exact knowledge of the spectrum. In what follows we shall associate our 
Fi(r, t) with the free energy F 1 (r,T), since its temperature dependence is 
not so steep as that of Ui(r,T) in this region and this discussion will be of 
qualitative character, leading detailed discussion of the spectrum to future 
publications. 
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3 Properties of D\(x) and F±(r,T) 



The correlator D\(x) was measured on the lattice jTH] both below and above 
T c , and decays exponentially with M 1 « 1 4- 1.5 GeV (in the quenched case). 
At the same time Di(x) can be connected to the gluelump Green's function, 
and the corresponding Mi for the electric correlator Df(x) is Mi ~ 1.5 GeV 
at zero T j2H]- Moreover in a recent paper [IQ Df(x) was found analytically 
for T < T c , and can be represented symbolically as a sum, with perturbative 
part acting at small x, 



Di(x) 



Df ert \x) + D?< <(.<) 



TTX H 



(14) 



and the nonperturbative part having the asymptotic form 



\x\ 



Mi\x\ 



+ 0(a% Ai = 2C 2 a s a adj Mi, x > I/Ml (15) 



As will be argued below, the form of D\(x) (|14j) does not change for T > T c , 
however the mass Mi and A± may be there different. 

Using the asymptotics ()15|) in the whole x region for a qualitative esti- 
mate, one has 



TO 



(np) 



(r,T)=AiJ^(l-uT)duJ ^^-^ 



Mi JO 



1/T, 



vT)dv[e 



-uMi 



^7-2+1/2 Mil 



M 2 



Ki(Mir)M ir - -^e 



e v, ,-„ AJ (16) 
Mir (l + Mir) +0(e" Ml / T )" 



Finally the Polyakov loop exponent is 



exp 



n (np) (o oi 
2T 



Vi np \oc) 



Ai_ 

M\ 



T 

Wi 



:i- e 



-Mt/Ts 



(17) 



One can see from (HID that u/ np) (oo) is finite and is of the order of few hun- 
dred MeV in the interval < T < 1.5T C . At small r from (JTfij) Vi m ' " 
const ■ r 2 . The total Vi 

bution at small r, which to the order 0(a s ) is 



y( n p) _|_ ^(p er *)^ cori ^ a i ns a i so perturbative contri- 



7T 



r/vfl-I/T) 



z/ 2 + 



y i (pert) (oo)+\/ 1 (C ' ) (r,T) 
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V} C \r,T) = —^±f(r,T), f(r,T) = 1 - - arctan(rT) - — ln[l + (rT)~ 2 ] 
r 7i 7i 

(18) 

From (fTHj) it is clear that Vi Pert \oo) is divergent and should be renormal- 
ized, Vi (oo) ~ 2C ^ Qs f~ — Tina) , a — > 0. Since the dominant divergent 

part is T-independent and (r) ~ r 2 at small r, one can renormalize 
matching Vi(r, T) with the Coulomb interaction at small r, as it was done in 
0HD for F;(r,T). 

As a result in the renormalized Vi Pert \r,T) the term V^ ert \oo) can be 
put equal to zero, and we shall use it in what follows. 

At this point we are able to compare Vi(r,T) with the lattice data for 
F\(r,T) at T > T c . In Fig.l we compare the lattice data for Fi(r,T) taken 
from [12] for T = 1.05T C , 1.2T C and 1.5T C with the potential V\{r, T) in the 
form (fTT)|) parametrizing Mi and a(T) = -pj- in it as 

T -T c 

Mi = const, a(T) = a — c — — — -, a s = 0.3 (19) 

J- c 

and find that Mi = 0.69 GeV and a = a{conf) = 2C 2 {f)a s a adj S 0.432 
GeV 2 , c = 0.36 provides a good agreement with the data points at 1.5T C > 
T > T c , while a(T) in (|19|) smoothly matches at T = T c the amplitude of 
the gluelump Green's function JH]- One can see that the behaviour of the 
total Vi(r, T) = V^(r, T) + v} np \r, T) which has a Coulomb part at smaller r 
and saturates at Vi = Vj np \oo) is qualitatively very similar to the behaviour 
of Fi(r,T) as a function of r. We also compare in Fig.2 our results with 
lattice data JT] for the Polyakov loop (JT7j) and find reasonable agreement. 
It is clear that both Fig.l and Fig.2 are qualitative illustrations, and for 
quantitative comparison one needs knowledge of excitation spectrum and 
analytic or lattice predictions for a(T), Mi(T) which will be given elsewhere 
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Figure 1: A comparison of behavior of V x (r,T) = v} np \r,T) + v}°\r,T) 
Eqs.(16),(18),(19)(solid lines, T/T c = 1.05; 1.2; 1.5 from above), with the 
singlet free energy Fi(r,T) measured in ref. ^2] (filled circles.) 
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Figure 2: The Polyakov loop exponent Vj(oo, T) as a function of T (GeV) 
from Eqs.(17),(19)(solid line), in comparison with lattice data for Fi(oo,T) 
in ref. [TT] (filled circles). 

4 Bound states of qq in the deconfinement re- 
gion 

In the recent lattice studies sharp peaks have been found in the spectral func- 
tion of cc system |HJ- [XO] , which can be associated with the quark- ant iquark 
bound states surviving at T > T c . 

To understand qualitatively whether the interaction V\ (r, T) can support 
bound states, one can use the Bargmann condition [27| for monotonic attrac- 
tive potentials 



where 2m = m c , and U(r, T) = Vi(r, T)— Vi(oo, T) which yields the condition 
for the bound ^-states. Taking m c = 1.4 GeV, Mi = 0.6 GeV, one can deduce 
that v} np \r,T) can support bound states in some interval of temperature 
T c < T < Td where To ~ 1.5 -4- 2T C and exact value depends on terms of 




(20) 
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the order 0(e~ Ml ^ T ) and therefore is beyond the scope of the present paper. 
This conclusion roughly agrees with lattice calculations in [H]-jH] and with 
the calculation done in [2E] and [2~9"j . 

One can consider gluons and the gg system in the same way as it was 
done for the cc system. To this end one should first multiply Vi(r, T) found 
earlier in (JHCHD by = §, thus defining V^\r,T) = \V X (r,T). This 

potential can be considered as the interaction kernel in the Hamiltonian of the 
gg system as it was done in [20]. This Hamiltonian with the introduction of 
the einbein gluon mass fi g = (y/p 1 ) has the form of the Schroedinger equation 
and the condition ([27?)) can be applied. Since /i 9 ~ 0.6 GeV j3U], an additional 
(with respect to quarks) factor for the nonperturbative part in (j2*U|) is ~ 
0.96. Hence two-gluon glueballs should be formed in approximately the same 
temperature interval as charmonium states. In lattice calculations [31] scalar 
glueball was studied below and around T c , and its width is increasing with 
T. Now we come to the baryon case and using the general formalism [32] to 
represent the triple Wilson loop of trajectories of 3 quarks and of the string 
junction in terms of field correlators D and D\. From [22] one has for 3 static 
quarks at distances ri, from the string junction 

V^(r 1; r 2) r 3 ; T) = £ V D (r t , T) + W - r v T) (21) 

i=l i>j 

where V\{r, T) is given in (|T6*j) . ()18J) . In the deconfinement phase when = 
R, i = 1, 2, 3, one obtains Vz q = |Vi(v / 3-R)- For the perturbative part one has 
from (J2TD 4 C) = lE l > J V 1 (C \r lJ ,T)\ rij=R . 

For the nonperturbative part from Eq. (|21j) it follows that V^g P \R — > 
^V} np \oo). One can check that this prediction and the general form 
of V 3q = ^Vi(y/3R) as function of R is supported by the recent measurement 
of singlet free energy of the 3Q system in j2H] at T > T c . Thus it is of interest 
to measure the spectral functions of baryons at T > T c in the same way as 
it was done for mesons. 

5 Summary and conclusions 

Citing the 1991 paper |5] when the magnitude of Di was not exactly known 
"...Using an exponential parametrization for Di, we can find Df with param- 
eter values which satisfy the condition for the appearance of levels. In this 
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case e(r) (our Vi(r, T)) is a well with a behaviour e(r) ~ r 2 as r — > and 
e(r) — > const > as r — > oo. The quark and antiquark are thus bound but 
there exists a threshold e(oo) above which quarks fly apart, each acquiring a 
nonperturbative mass increment 8m = |e(oo)...". In the present paper this 
picture was further substantiated and quantified using lattice and analytic 
knowledge on D\. Comparing to recent lattice data in J3| it was shown 
that this picture is qualitatively supported by data, and new proposals have 
been done for searching the glueball and baryon systems at T > T c . 

The results obtained on the lattice 0-^21 and in the present approach es- 
tablish a new picture of the QCD thermodynamics at 1.5T C > T >T C widely 
discussed in [2E]- As a new feature compared to the works [2E] the main 
emphasis in this paper is done on the selfenergies (|Vi(oo) = e q for quarks 
and |Vi(oo) = e g for gluons) which are large (Vi(oo,T c ) > Fi(oo,T c ) = 600 
MeV for rif = 2 [S3]) and cancel each other at small distances for white 
bound states, like qq, (gg)i, (qgq)i, (qg---gq)i etc. In contrast to that colored 
states are higher in potential and mass by several units of e q and e g and 
are suppressed by the corresponding Boltzmann factors. As a result in this 
region white bound states of quarks and gluons are energetically preferable, 
while individual quarks and gluons acquire selfenergies, so that the thermo- 
dynamics of the system resembles that of the neutral gas, and for higher 
temperature T > 1.5 4- 2T C a smooth transitoon to the "ionised" plasma of 
colored quarks and qluons possibly occurs. 

This new state of the quark-gluon matter should be taken into account 
when considering ion-ion collisions. For more discussion of the thermody- 
namics above T c see j2H] and refs. therein. 

It was noted before [3I| that the behaviour of the free and internal energies 
above T c , with a bump around T ~ 1.1 -j- 1.2T C in e ~^ p can be explained 
if gluons are supplied with the nonperturbative mass term of the order of 
0.6 GeV, while for higher T this mass is less important. This can be easily 
understood now taking into account the value of |Vi(oo, T) and its decreasing 
with growing T. In this way the nonperturbative dynamics in the form of 
correlator D\ can explain the observed dynamics of the deconfined QCD. 

A more detailed analysis of bound states requires explicit calculation of 
QQ and 3Q bound states taking into account spin splitting in the mass of 
P-wave charmonia and quasi-degeneration of spectra of light qq V, A, S, PS 
states observed in [T3j . Here spin-dependent forces are different from the 
confining case, since only the correlator D\ contributes, and one can list the 



12 



corresponding terms in [T71 EJB]. It is interesting to note, that due to the 
vector character of V\{r,T), not violating chiral symmetry, bound states of 
massless quarks should exhibit parity doubling. All this analysis is now in 
progress [3*7] . 
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Appendix 1 

Derivation of the Polyakov loop correlator 

We give here two different derivations of Fi(r, T), i — 1,8. The first is 
based on the correlator of two concentric Wilson loops, derived in [20], in 
which case Fi$ are expressed in terms of surface integrals of field correlators 
F>v, v , p \{u,v) 




(A1.1) 



In this way one obtains for two oppositely directed Polyakov loops from 
Eqs. (29), (23), (20) of [H| 

F 1 (r,T)/T=~I(S 12 ,S 12 ) (A1.2) 

F 8 (r, T)/T = ±I(Si, St) + h(S 2 , S 2 ) + jJ-^Ifa, S 2 ) (A1.3) 

Here S\ is the surface on the cylinder with circumference 1/T extending 
from the loop 1 at coordinate x in the direction y to infinity, the surface 5*2 
is also infinite surface from the loop 2 at coordinate y in the same direction 
(the answer does not depend on the choice of this direction). 
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The surface Si 2 lies on the cylinder between the loops 1 and 2. Note that 
surface orientation in (|A1.1J) is fixed to be same. Calculation of F\ according 
to (jA1.2j) reduces to that of the Wilson loop and yields 

F 1 (r,T)=V 1 (r,T) + V D (r,T) (MA) 

where V\ (r, T) is given in (jlj) and Vo is 

rl/T 

du{l - vT) 

lo 



V D (r, T) = 2 du(l - uT) f (r - 0<%D(y/e + z/ 2 ) (A1.5) 
Jo Jo v 



Calculation of F$ is more subtle. To this end one can use connection of 
D\(x) to the gluelump Green's function G^ v (x) = 5^ U N C (N 2 — l)f(x 2 ) [H] 
D x {x) = - 2 -f c (NZ - 1)^1, f(x 2 - 0) ~ ^ and inserting this into (4), 
one has 

2(jy-2 _ i\ r i/ T 

Vi (r, T) = 2_L^ L / dv (i^ T )[f{i?)-f{r 2 +p 2 )] = V\(oo, T)+v ex (r , T) 

iv c Jo 

(A1.6) 

One can see that Vi(oo,T) = Vq + Vq is the sum of equal selfenergy 
parts of Q and Q, while w ex (r, T) describes interaction due to one gluelump 
exchange between Q and Q. 

Note that Vi(0, T) = 14(00, T)+ v ex (0,T ) = and w ex (oo, T) = 0. There- 
fore v ex appears only in I(Si, S 2 ) in ()A1.3|) and one should restore there the 
original (opposite) orientation of loops L x and Ly to get the correct sign of 
v ex (the same sign and factor appears in the second derivation below). 

From (jAl.l|) one obtains for I (Si, Si) 

l -I(S 2 , S 2 ) = V D (r*) + V Q , h(S l , St) = \l(S 2 , S 2 ) + V d (t), 

I(S U S 2 ) = -v ex (r,T) + I(S 2 ,S 2 ) (A1.7) 
As a result one can write using (|A 1.3 J) for N c = 3 

F 8 (r,T) = \v D (r\T) + Vq + Vq — l -v ex (r,T) + 

+ V D (r,T) = ^1(00, T) - ^(r,T) + ^V D (r*,T) + V D (r,T) (A1.8) 

In (|A1.8|) the value of r* is infinitely large, when one neglects the valence 
gluon loops, as it is done everywhere above. In this case Vjj(r*, T) — > 00 and 
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the term exp(— F$/T) vanishes in the confinement region. This is in line with 
the strong damping of the adjoint Polyakov loop in this region observed on 
the lattice and with the persistence of Casimir scaling for adjoint static 
potential for < r < 1.2 fm found on the lattice (see discussion and refs.in 



The correction due to the gluon determinant, producing additional gluon 
loops becomes important for r > 1.2 fm (see discussion in the second ref. in 
and makes finite the value of Vo{r*,T) « ar*,r* « 1.2 fm/2 =0.6 fm. 



exp(-F 8 (r,T)/T) = L adj (T)exp 



V D (r,T) - IV^TY 
T 



L adj (T) = exp(-^y 1 (oo,T) - ~V D (r*,T)). (A1.9) 

In (jA1.9|) the effects of loop-loop interaction and of the total (adjoint) 
loop are separated. Physically the result (jA1.9|) can be easily understood: in 
absence of the internal interaction one has do with the adjoint Polyakov loop, 
which strongly changes around T c , namely L a( ij{T < T c ) is much smaller than 
L a dj(T > T c ). 

The behaviour similar to (|A1.9|) was observed on the lattice, see Fig. 3 of 
second ref. jTUj and ref. jMj. Here one observes linear growth in r of F 8 (r, T) 
below T c , and repulsive Coulomb behaviour from |Vi(r, T) = 

In the second derivation one is connecting two loops using parallel trans- 
porters and using the completeness relation 

(WifWa = J-j-ficafatipKxz + 2tp 2Ctl tp ia2 (ALIO) 

so that the second term produces the adjoint Wilson loop on the cylinder 
surface. 

exp(-F 8 /T) = 2ir((£/(x,X; 0)t a U(X, y, O)L(y) x 

x U(y, X, t)t a U(X, x, t)L + (x)) (Al.ll) 

This can be compared to the approach, suggested in j^HJ. Our results 
however are different from those of [SH] m that both perturbative and non- 
perturbative interactions in F\ and F§ are different (and calculable through 
Di). Now (jAl.llJ) can be rewritten using cluster expansion and nonabelian 
Stokes theorem, which finally results in the same equations as in (|A1.9|) . 
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Alternatively one can use the technic exploited in [3D] to separate the contri- 
butions of perturbative exchanges from the nonperturbative confining terms. 

For the first ones one commutes as in [SH] the color generators t c of ex- 
changed gluon (gluelump) with t a flAl.ll|) acording to the equality t c t a t c = 
—t a /2N c which finally gives in (jAl.lljl the adjoint Coulomb interaction ^ — > 
— |Vi(r), while the selfenergy parts and confining terms arise from sequences 
t c t c 't a — > 5 cc /t a and do not change sign. In this way one arrives at the same 
answer as given in (jA1.9j) . 
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